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ABSTRACT 

We present high-resolution spectroscopy and images of a photodissociation region (PDR) in M16 
obtained during commissioning of NIRSPEC on the Keck II telescope. PDRs play a significant role in 
regulating star formation, and M16 offers the opportunity to examine the physical processes of a PDR 
in detail. We simultaneously observe both the molecular and ionized phases of the PDR and resolve the 
spatial and kinematic differences between them. The most prominent regions of the PDR are viewed 
edge-on. Fluorescent emission from nearby stars is the primary excitation source, although collisions also 
preferentially populate the lowest vibrational levels of H2. Variations in density-sensitive emission line 
ratios demonstrate that the molecular cloud is clumpy, with an average density « = 3x 10 5 cm~ 3 . We 
measure the kinetic temperature of the molecular region directly and find Tjj 2 = 930 K. The observed 
density, temperature, and UV flux imply a photoelectric heating efficiency of 4%. In the ionized region, 
rii = 5 x 10 3 cm -3 and Tjjn = 9500 K. In the brightest regions of the PDR, the recombination line widths 
include a non-thermal component, which we attribute to viewing geometry. 

Subject headings: infrared: ISM: lines and bands — ISM: individual (M16) — molecular processes 



OS 

(N 
O 
O 

o 
6 



1. INTRODUCTION 

Star formation occurs in dense regions of the interstel- 
lar medium, and ao a rooult, the onvironmonto of early 



cal emission lines ( Hester ct al. 1996| ). The young cluster 
NGC 6611 contains many massive stars, which illu minate 
the column s with total FUV flux 16 ergs^ 1 cm' - 2 (lATIeil 



typo at Bra ofton contain molooular oloudo i In addition 



to creating an H II region, where hydrogen is predomi - 



1996) 



et al. 1999) from a distance of about 2 pc (Hester et al 



We simultaneously observe the spatial and spectral re- 
lationship of the excited molecular and ionized phases of 
the PDR. Based on these observations (described in §2), 
we determine physical conditions in §3, and summarize the 
conclusions in §4. 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained images and high resolution longslit spec- 
tra of the northernmost column of M16 using NIRSPEC, 
the facility n ear-infrared spectro meter on the Keck II 10- 
m telescope ( McLean et al. 1998| ). Figure |l|a identifies the 
observe d region on the Hu bble Space Telescope Ha image 
of M16 ( Hester et al. 1996 ) with contours of integrated CO 
intensity ( pound 199S ) ■ We used the slit- viewing camera 
to image M16 on 1999 July 6 (UT), obtaining 9 exposures 
of 60 s each. In addition to the H2 1-0S(1) filter, (bandpass 
2.110-2.129 fim) the NIRSPEC-6 filter (1.558-2.315 //m) 
was in the optical path. We combined several observations 
of a blank field in the NIRSPEC-6 filter alone to use as 
a flat field and scaled this for sky subtraction from this 
crowded field. Figure shows the resultant 2.12 /im im- 
age and the spectroscopy slit. 

We employed the NIRSPEC-6 blocking filter and se- 
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nantly ionized, stellar ultraviolet radiation creates a pho 
todissociation region (PDR) when it interacts with molec- 
ular material. The feedback between extant stars and the 
nearby clouds in PDRs is important in regulating subse- 
quent episodes of star formation. The UV flux determines 
the ionization fraction in molecular clouds, which in turn 
sets the ambipolar diffusi on rate and thereby the rate of 
low- mass star formation ( |McKee 198S ). 

Hollenbach & Tielens (1997, 1999) comprehensively re- 
view the physical and chemical processes of PDRs. Far- 
ultraviolet (FUV) photons (6 < hv < 13.6 eV) dissociate 
molecules and excite the gas on the surface (Ay < 3) of 
the molecular cloud. The result is spatial stratification of 
molecular, atomic, and ionic components of a given ele- 
ment. At higher densities (n> 10 4 cm~ 3 ), collisional exci- 
tation also populates the excited states ( [Sternberg fc Dal 
garno 1989|) 



The columns of M16 (the Eagle Nebula) are PDRs. 
Most of the mass in the nebula is in molecular clouds, 
and a variety of molecular transiti ons are observed directly 
( |Pound 199^ ; |White ct al. 1999|) H 2 is observed on the 
cloud surfaces ( Allen et al. L99S| ) , and at the interface 
with the H II region, the columns are prominent in opti- 
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Fig. 1. — (a) The HST He* image of M16 (Hester ct al. 199C), with contours of integrated CO intensity {solid lines; Pound 199£) and the 
location of the NIRSPEC slit-viewing camera 2.12 fim image (dotted line) overlaid. The Ho image is scaled with a square-root stretch, and 
the CO contours are scaled linearly, from 50 to 300 Jykms . (b) Northern column of M16 at 2.12 /im (including H2 1-0S(1)), observed 
with the NIRSPEC slit- viewing camera. Right ascension and declination are marked with respect to the observed center of the spectroscopy 
slit, which is overlaid. The UV flux comes from the northwest, which corresponds to the top of the slit in subsequent figures. The excited 
molecular emission is a surface phenomenon, arising only at the edge of the molecular cloud. 



lected echelle and cross disperser positions to detect sev- 
eral H2 rovibrational transitions, Br7, and He I A2.06 fim 
in our spectroscopic observations. We used the (X'432 x 24" 
slit at position angle 135° for resolution A/AA = 26,000, 
where AA is the observed FWHM of an unresolved emis- 
sion line. We obtained six exposures of 300 s each on the 
night of 1999 April 30 (UT), nodding along the slit be- 
tween each exposure. 

We subtracted a median dark frame from each image, 
divided by a flat field, and interpolated over deviant pix- 
els, removing both cosmic rays and bad detector pixels. 
The spectra required rectification in both the spectral and 
spatial dimensions. We used the OH night sky lines and 
wavenumbers tabulated by Abrams et al. (1994) for wave- 
length calibration and the continuum emission of the tel- 
luric standard for spatial rectification. We observed HD 
f6f056 and reduced the spectra similarly, then modelled 
this B1.5V star as a T = 24,000 K blackbody and inter- 
polated over the stellar Br7 feature. We divided the M16 
data by this spectrum to correct for atmospheric absorp- 
tion and normalized the stellar spectrum from this K = 5.4 
standard at 2.2 /zm to flux calibrate the data. 

Figure || contains the resultant two-dimensional spec- 
tra of several bright lines. We have subtracted off-source 
emission from each of these spectra. For the molecular 
emission, the background is near the edge of the slit to- 
ward the ionizing source, while for Br7, the background 
lies away from the ionizing source. We list the complete 
set of detected lines, their vacuum wavelengths, their total 
fluxes in the slit, and their peak fluxes in Table [l] 



3. PHYSICAL CONDITIONS 

3.1. Geometry 

The gross structure of the photodissociation region is 
evident in the images (Fig. and two-dimensional spec- 
tra (Fig. ||). Fewer stars are observed directly through 
the cloud than in the surrounding region (Fig. |I}d), im- 



plying Ak > 1- From the center of the molecular cloud 
out to the photoionizing source, distinct regions of qui- 
escent molecular, excited molecular, atomic, and ionized 
material appear. The FUV photons do not fully pene- 
trate the molecular cloud, so the excited H 2 resides only 
on the cloud's surface. The optically-thin ionized emis- 
sion is brightest in the regions of highest density, near the 
cloud, and it is spatially more extended than the H2 emis- 
sion. The slit is aligned toward the ionizing source and 
covers several edges on the irregular surface of the cloud. 
Thus, we observe several distinct regions of photoevapo- 
rative flow. The greatest intensities are detected in the 
regions that are viewed edge-on. 

In these observations, each of the bright regions within 
the PDR consists of spatially segregated molecular and 
recombination line emission. The two brightest regions at 
slit positions — 1'/2 and — 3'.'2 in the 1-0S(1) emission, for 
example, are physically associated with the ionized emis- 
sion at 7'.'1 and 3'.'6, respectively (Fig. ||d). (Note, how- 
ever, that the more spatially extended Br7 peak at the top 
of slit includes some contribution from a nearby surface 
north of the slit, which we do not observe in H2.) Mod- 
erate emission peaks appear at 2'.' 8 and 1'.'8, and fainter 
surfaces at -9'.'3, -7'.'9, and 5'.'3 in the 1-0S(1) line. The 
first two of these are blended in the spatially extended 
Br7 emission, while the latter correspond to local emission 
maxima at — 2'.'8, — 0'.'2, and 11'.' 3, so the projected sepa- 
rations of the H2 and Br7 emission regions in the plane 
of the sky range fro m 1.8 to 2.5 x 10 17 cm at the 2000 



pc di stance of M16 (Humphreys 1978; Hillenbrand et al 



1993). The transition from t he atomic to molecular zon e 



typically occurs at Ay ~ 2 (Ticlens fc Hollcnbach 1985), 



or in terms of column density, Nh 
the average gas density in the bright regions of the PDR 
is 1.6-2.2 x 10 4 cm" 3 . 

We determine the density in the ionized region and use 
the molecular line width (discussed below) to calculate 
the average density in the molecular region of the PDR. 
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Bertoldi & Draine (1996) show m = Sl v /(47r R 2 qci), where 
Si, y is the Lyman continuum flux, R is the distance to the 
ionizing source, q accounts for attenuation of the ioniza- 
tion in the evaporative flow, and c~ is the isothermal sound 
speed in the ionized mediu m. We determ ine the value of 
q ~ 4.0 x l(r 14 ySl~F/i? flBcrtoldi 19891 ), wh ere r is the 
cloud radius. We find the density of the ionized region 
5.2 x 10 3 cm -3 , for g = 12 kms -1 (calculated for 



T Hu 4= 9500 K, below) and r = 25". For pressure bal- 
ance across the ionization front, n ~ 2n,icf/(Av/2.35)' ! , 
where Av is the FWHM of the molecular line width, so 
in the molecular region, n ~ 3.2 x 10 5 cm -3 . This is signifi- 
cantly greater than the density calculated above, which in- 
tegrated over contributions from the ionized to the molec- 
ular region. Here, n is determined in the molecular zone 
alone, which should be denser. 
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Fig. 2. — Two-dimensional spectra of bright lines, (a) H2 1-0S(1 
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fim, scaled linearly from to 1.8 X 
2 . (c) Br7, 2.17 (mi, scaled linearly from 



(b)H ? 1-0S(0), 2.22 
10 -15 ergs -1 cm™ 2 

to 3.6 X 10 -15 ergs -1 cm -2 arcsec -2 . (d) Contours of Br7 (solid 
lines) and H2 1-0S(1) (dotted lines) emission, in logarithmic intervals 
from 0.6 to 3.4 X 10 -15 ergs -1 cm -2 arcsec -2 . The Br7 emission is 
broader and closer to the ionizing source (located at the top of the 
slit) than the molecular lines. The two brightest regions in Br7 at 
slit positions 3'.' 6 and 7'.'1 correspond to the H2 peaks at —3'.' 2 and 
— 1"2, respectively. 



3.2. Kinematics 

The various gas phases exhibit distinct kinematic struc- 
ture. We measure the line profiles at slit position 3'.'2 in 
1-0S(1) emission, which is clearly distinguished in the Br7 
and He emission at 3'.'6. Because we observe this PDR 
nearly edge-on, the central velocities of the H2 and ion- 
ized material are not significantly different; the average 
motion occurs in the plane of the sky. We measure vlsr — 
29 kms" 1 for this 1-0S(1) feature and vlsr = 28 kms" 1 in 
Br7. InH 2 , Aw = 4.6±0.1 kms" 1 , Av = 29.0±0.2 kms" 1 
in Br7, and Av = 18.7 ±0.2 kms" 1 in He I A2.06^m, after 
correcting for the instrumental profile. In the molecular 
gas, the velocity width corresponds to a kinetic tempera- 
ture Tjj 2 — 930 ± 50 K. This value is reasonable near the 
edge of the cloud in the H ? zone where maximum tem per- 
ature of the PDR occurs (Fielens & Hollcnbach 1985). 



Accounting for the mass difference between H and He 
demonstrates that these line widths are not purely ther- 
mal, however. Assuming that they each consist of a 
macroscopic component that is independent of molecu- 
lar weight in addition to the thermal component, we find 
Avmacro = 12.2kms -1 . For Br/y, Av t herm = 25.6kms" 1 , 



and T HlI = 14000 ± 230 K. This temperature is hotter 
than typically observed in H II regions, but several factors 
tend to increase Tjjn at the PDR-H II interface: increased 
photoelectric heating near a greater population of grains; 
the harder radiation field, which reduces the contribution 
of significant low-ionization coolants, such as + ; and in- 
creased density, which reduces the radiative cooling rate by 
incre asing the collisional de-excitation rate (e.g., 3pitzei 



1978) 



The non-thermal recombination line widths may also 
be attributable to bulk motion of material flowing off the 
cloud surface. In the bright region discussed above, we 
view the PDR nearly edge-on, along a line of sight that 
includes both rcdshiftcd and blucshiftcd velocity compo- 
nents of the material that evaporates off the curved surface 
of the cloud. When the view of the PDR is closer to face- 
on, the intensity is less enhanced, and only the blueshifted 
velocity component arises. This is the case for the minor 
emission peaks at the bottom of the slit, which occur on 
the near surface of the cloud. These central velocities are 
slightly lower (by approximately 7 kms" 1 ) than those of 
the brightest peaks, and the observed widths are narrower. 
At slit position — 2'.'8, for example, we find Av = 20.9 ±0.2 
and 10.1 ±0.2 kms" 1 in Br7 and He I, respectively. These 
both correspond to kinetic temperature Tjj n = 9500 ± 180 
K. If this is a more accurate measurement of the temper- 
ature at the edge of the H II region, then we have un- 
derestimated the contribution of the non-thermal velocity 
component above. While Tjjn at the interface is likely 
higher than the typical temperatures in the diffuse central 
portions of H II regions, as argued above, it may be less 
than 14,000 K. 

3.3. H% Excitation 

In a PDR, UV pumping is the primary source of H2 exci- 
tation. In M16, such fluorescent emission obviously occurs. 
A variety of lines that originate in high vibrational levels 
v are observed, including 7-50(4) and 3-2S(5) (Table 0), 
whereas thermal population of these high-u levels would 
dissociate the molecules. We detect emission from all of 
the strong transitions in our bandpass that Black and van 
Dishoeck (1987) predict, except those originating in levels 
v > 9 or lie in regions of poor atmospheric transparency. 

Collisional excitation is also important, however. The 
ratio of line intensities 7(l-0S(l))/7(2-lS(l))= 2-5, and 
the ave rage ratio over the slit is 4. As we found above 
in § 3.1 , this req uires n > 10 4 cm" 3 (Burton, Hollenbach 



fc Tielens 1990] ), so collisions aid in populating the v = 1 
level. We do not, however, observe line ratios as high as 
Allen et al. (1999), who find 7(l-0S(l))/7(2-lS(l))> 10 in 
similar regions with Fabry-Perot imaging. We expect our 
results to be more accurate, having better flux calibration 
and simultaneous sky measurements. 

Line ratios also directly measure excitation tempera- 
tures. Across the slit, we find from the 1-0S(1) and 2-lS(l) 
lines vibrational temperature T v n, = 3000-3800 K, with a 
median T vib = 3400 K. Using 7(1-0S(1))/7(1-0S(3)), we 
find rotational temperature T rot = 1600-2100 K and me- 
dian T rot = 1800 K. Comparing 7(1-0S(1))/7(1-0S(0)), 
T rot = 450-870 K, with a median T rot = 580 K. The 
relatively low value of T v n, and relatively high value of 
Trot (particularly in the brightest regions) distinguish the 
M16 PDR from purely fluorescent emission in a low den- 
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sity medium, where T v n, ~ 5500 and T rot ~ 1400 from 
J(1-0S(1))/I(1-0S(3)) are typical. 

We compare the line intensities of the three brightest 
regions with the models of Draine & Bertoldi (1996). Con- 
sidering all observed lines, all three regions are best fit with 
x/n — 0.1 cm 3 , and To = 1000 K, where x measures the in- 
cident UV flux (x = 1 is the Habing [1968] flux), and T is 
the temperature at the edge of the PDR. The observed line 
ratios clearly rule out models of a weak field (x < 10 3 ) and 
low density (n < 10 4 cm -3 ). Both x = 10 4 and x = 10 5 
are acceptable. For consistency with other density mea- 
surements and the UV field expected from members of 



NGC 6611, x = 10 4 and 



10 cm 3 is the preferred 



model. Although the grid of models is coarse, we iden- 
tify trends by contrasting these regions. For example, the 
variation of density-sensitive line ratios along the spatial 
extent of the slit demonstrates that the cloud density is 
not constant. The decreased ratios J(1-0S(1))/I(1-0S(0)) 
and I(1-0S(1))/J(2-1S(1)) at slit position -3'.'2 relative to 
slit position — 1'/2 show that the former region has lower 
density. 

In thermal equilibrium, the ratio of ortho-hydrogen 
(having odd rotational level, J) to para-H2 (having even 
J) is 3, and only collisions can change the ortho:para ratio. 
This abundance ratio is therefore a significant parameter 
in PDRs and affects the emission line ratios. We calculate 
the ratio of ortho-H2 and para-H2 from the 1-0S(1) and 1- 
0S(0) lines, finding ratios of 1.3-2.2, with an median value 
of 1.7. Because UV pumping of the damping wings of op- 
tically thick lines populate the excited levels, the expected 
ratio of column densities in the vibrationally exc ited levels 
is 1.7, while the true ortho:para abundance is 3 (Sternberg 
& Ne jfeld 1999| ). 

The photoelectric effect on interstellar grains is likely to 
be the primary heating source in PDRs, but its exact pro- 
cess is uncertain. We determine the efficiency of convert- 
ing FUV flux to photoelectric heating, e, assuming ther- 
mal equilibrium. Emission in the [O I] 63 /im and [C II] 
158 /mi fine structure lines cools the gas. For the mea- 
sured temperature and density, we calculate a net cooling 
rate A = 4.2 x 10" 17 erg s _1 cm" 3 in M16, with 96% due to 
[O I]. In the parameterization of Bakes & Tielens (1994), 
the heating function T = 10~ 2 ne^ergs" 1 cm" 3 . Thus, 



we find e = 0.042, for X = 9700 (lAllen et al. 1999]). Th e 
theoretical prediction, e = 0.016 ( [Bakes fc Tielens 1994 ), 
is a function of XvT/n e , where n e is the electron density. 
We have assumed that the ionization of carbon provides 
the electrons; increasing n e would increase e, without sig- 
nificantly altering the cooling rate. Alternatively, decreas- 
ing the gas to dust ratio would allow a lower efficiency in 
equilibrium with the cooling rate we predict. 

4. CONCLUSIONS 

Spectroscopy of M16 with NIRSPEC offers the advan- 
tage of simultaneous observation of both the molecular 
and ionized phases of its photodissociation regions. At 
high resolution, we detect the spatial and kinematic sig- 
natures of photoexcitation of molecular material. Most of 
the emission is due to distinct regions of photodissociation 
on the molecular cloud's surface that we observe edge-on. 
The spectrum includes several lines that originate in high 
vibrational levels, which demonstrates that fluorescent ex- 
citation predominates. The density n ~ 10 5 cm" 3 and 
varies among the emitting regions, so collisions also pref- 
erentially populate low- lying vibrational levels. We mea- 
sure kinetic temperature Th 2 — 930 K in the molecular 
gas, which requires efficient photoelectric heating near the 
cloud's surface. In the spatially-segregated ionized region, 
Tmi = 9500 K, and the velocity widths of the brightest 
regions include a macroscopic contribution. 
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Table 1 
Observed Lines 



Transition Wavelength 


Total Flux 


Peak Surface Brightness 


[Mm] 


[ergs _1 cm~ 2 ] 


[erg s ~ 1 cm ~ 2 ar csec ~ 2 ] 



H 2 



7-50(4) 


1.94343 


3.9 ±0.2 x 10~ 16 


1.9 ±0.4 x 


10" 


■16 


2-lS(5) 


1.94486 


1.1 ± 0.01 x 10~ 14 


3.2 ±0.2 x 


10" 


■15 


l-0S(3) a 


1.95756 


5.6 ± 0.02 x 10~ 14 


4.8 ±0.2 x 


10" 


■15 


3-2S(5) 


2.06556 


4.5 ± 0.06 x 10~ 15 


3.7 ±0.6 x 


10" 


■16 


8-60(4) 


2.12155 


7.4 ±0.2 x 10~ 16 


1.2 ±0.3 x 


10" 


16 


1-0S(1) 


2.12183 


6.2 ± 0.02 x 10~ 14 


5.3 ±0.2 x 


10" 


■15 


3-2S(4) 


2.12799 


1.9 ±0.04 x 10" 15 


2.4 ±0.5 x 


10" 


16 


1-0S(0) 


2.22329 


2.4 ±0.01 x 10~ 14 


1.6 ±0.1 x 


io- 


■15 


2-1S 1 


2.24772 


1.6 ±0.01 x 10" 14 


1.3 ±0.1 x 


io- 


■15 


Recombination Lines 










Br7 


2.16609 


1.3 ±0.003 x lO" 13 


3.7 ±0.2 x 


io- 


■15 


Br5 a ' b 


1.94509 


2.3+l Q2 x 10- 14 


1.3 ±0.1 x 


io- 


■15 


He I 


2.05869 


9.7 ±0.03 x 10" 14 


3.9 ±0.2 x 


io- 


■15 


He I 


2.11259 


2.6 ±0.05 x 10~ 15 


1.9 ±0.4 x 


io- 


■16 



a Correction for atmospheric transparency is large, around 2 
b Missing flux due to poor atmospheric transparency 



